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ABSTRACT 


It  has  been  experimentally  observed  for  some  time 
that  certain  tearing  modes  in  plasmas  may  be  suppressed 
if  the  plasma  rotates  in  a  preferred  direction.   In 
this  paper  we  treat  the  m=0,  finite  wavelength  tearing 
mode  in  cylindrical  geometry  for  a  reversed  field  plasma 
equilibrium  and  show  that  by  generalizing  Ohm's  law  to 
include  Hall  current  terms,  we  are  able  to  explain  this 
effect  of  rotation  on  tearing  modes.   Our  results  agree 
qualitatively  with  earlier  analysis  and  numerical 
simulations  of  A.  Kaleck.   We  also  show  that  our  results 
are  sensitive  to  the  position  of  the  outer  conducting 
wall,  and  for  wall  positions  sufficiently  close  to  the 
plasma  vacuum  interface,  tearing  modes  may  be  quenched 
when  the  rotation  reaches  a  critical  value.   These  results 
follow  from  a  boundary  layer  analysis  and  numerical  inte- 
gration of  the  boundary  layer  equations. 
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1.   INTRODUCTION 

The  purpose  of  the  present  work  is  to  re-condider  a  well-known 
resistive  MHD  instability  by  considering  the  effect  of  a  more  gen- 
eral form  of  Ohm's  Law.   Specifically,  we  treat  the  m=o,  finite  wave- 
length tearing  mode  in  cylindrical  geometry  for  a  reversed  field  plasma 
equilibrium   typical  of  those  which  are  produced  in  reverse-bias  theta 
pinch  experiments.   (BODIN,  1963;  KALEK  et  al.,  1968;  EBERHAGEN  and 
GROSSiMANN,  1971;  ES '  KOV  et  al.,  1978;  LINFORD  et  al.,  1978).   The 
motivation  for  the  present  study  stems  from  early  experiments  which 
showed  a  curious  and  unexplained  absence  of  an  expected  tearing  mode 
(EBERHAGEN  and  GROSSMANN,  1971).   Reverse  bias  theta  pinch  experi- 
ments generally  produce  a  closed  field  configuration  with  reverse 
trapped  magnetic  flux;  this  occurs  due  to  magnetic  line  reconnection 
at  the  cold  ends  of  the  pinch  coil  where  the  resistivity  is  high. 
The  formation  of  the  closed  region  is  also  associated  with  a  strong 
axial  contraction  as  the  plasma  attempts  to  find  an  equilibrium. 
The  result  is  an  elongated  (length  to  width  ratio  >  10)  reversed- 
field  plasma  with  a  bounding  closed  separatrix,  which  according  to 
the  simplest  classical  picture  should  be  further  susceptable  to  tear- 
ing modes,  m=o  and  k   ^  0 ,  which  would  tend  to  break  the  column  into 
a  number  of  cells  each  having  reversed  trapped  flux.   The  number  and 
size  of  the  cells  as  predicted  by  the  simple  linear  theory  varies 
from  experiment  to  experiment  since  the  linear  dispersion  relation 
depends  sensitively  on  plasma  profiles.   The  classical  tearing  mode 
theory  due  to  FURTH,  KILLEEN  and  ROSENBLUTH  (1963)  carried  out  original- 
ly in  slab  geometry  yields,  for  cylindrical  geometry  (KALEK,  1972) 
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the  prediction  that  most  of  the  experimentally  produced  reversed 
field  theta  pinch  plasmas  should  break  up  into  reversed  field 
cells.   Experimental  results  of  EBERHAGEN  and  GROSSMANN  (1971)  and 
more  recently  by  LINFORD  (1978)  show  that  the  expected  tearing  mode 
does  not  occur  and  the  plasmas  in  question  exist  for  some  tens 
of  linear  growth  times.   Instead  of  the  expected  tearing  a  pro- 
nounced rotation  of  the  plasma  column  is  observed;  this  rotation 
increases  with  time  until  the  rotational  frequency  reaches  a 
critical  value  at  which  time  the  plasma  goes  catastrophically 
unstable  to  an  m=2  mode.   In  an  attempt  to  explain  the  absence 
of  the  tearing  modes  KALEK  (1972)  simulated  the  plasma  behaviour 
with  a  linearized  initial  value  MHD  code  assuming  the  plasma 
column  to  be  infinitely  long  with  no  longitudinal  spatial  de- 
pendence.  Kalek's  model  originally  contained  the  resistive  MHD 
equations  of  motion  linearized  about  a  rigidly  rotating  plasma 
equilibrium;  results  from  the  numerical  code  showed  that  the 
inclusion  of  rotation  predicted  even  larger  growth  rates  than 
for  the  static  equilibrium  case.   The  addition  of  the  Hall  effect 
into  the  Ohm's  Law  yielded  results  which  were  not  inconsistent 
with  the  experimental  results,  namely  the  growth  rates  as  cal- 
culated were  at  least  an  order  of  magnitude  lower  than  what  the 
FKR  theory  would  predict;  this  statement  requires  qualification 
in  that  the  relevant  part  of  the  FKR  cylindrical  tearing  mode 
dispersion  relation  is  distorted  in  a  manner  to  suppress  the 
growth  rate.   The  long  wavelength  (ka  <_  1)  part  of  the  dispersion 
relation  is  affected  and  a  preferred  direction  of  rotation  is 
indicated  for  a  stabilizing  influence.   The  direction  of  the  ion 
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diamagnetic  drift  velocity  is  preferred  for  stability  and  ro- 
tation in  the  opposite  direction  results  in  enhanced  instability. 
Kalek's  numerical  results  lend  strong  support  to  a  beneficial 
interplay  between  the  effects  of  rotation  and  Hall  current  in 
the  suppression  of  tearing  mode  growth  rates.   The  calculation 
described  in  the  present  paper  attempts  to  analytically  in- 
vestigate the  m=o,  K   7^  0  tearing  mode  in  reversed  field  plasmas 
and  attempts  to  explain  the  numerical  results  of  KALEK.   Our 
results  are  consistent  with  Kalek's  and  we  find  a  stabilizing 
influence  of  rotation  combined  with  Hall  effect;  however,  we 
find  that  the  results  are  very  sensitive  to  profiles,  wall 
position  and  plasma  parameters  such  as  density,  temperature  and 
rotation  frequency.   The  specific  calculation  presented  in  the 
present  paper  deals  with  plasma  profiles  which  are  representative 
of  those  which  are  produced  in  reversed  field  theta  pinch  experi- 
ments. 

It  should  be  mentioned  that  the  use  of  Ohm's  Lav;  with  a 
Hall  term  has  also  been  used  to  model  the  effect  of  finite  par- 
ticle effects  in  treating  tearing  modes  since  the  Hall  term  is 
associated  with  a  coefficient  which  is  proportional  to  the  ion 
gyration  frequency  (COHEN  et  al.,  1977).   In  such  a  case  over- 
stable  modes  can  be  found.   In  the  problem  treated  here  the  m=o 
mode  does  not  admit  overstable  solutions;  this  is  not  explicitly 
pointed  out  in  the  text  but  follows  from  simple  energy  arguments. 

The  paper  is  organized  into  9  Sections.   The  governing  equa- 
tions are  presented  in  Section  2.   The  governing  equations  are 
non-dimensionalized  in  Section  2.   Solutions  of  the  governing 
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equations  in  the  resistive  boundary  layer  are  discussed  in  Sec- 
tion 3.  Solutions  outside  the  resistive  layer  are  discussed  in 
Section  4  and  matching  across  the  boundary  layer  is  discussed 
in  Section  5.  The  tearing  mode  ordering  with  a  small  parameter 
is  discussed  in  Section  6.  Sections  7  and  8  discuss  the  com- 
putation of  growth  rates  for  specific  equilibria.  A  final  sec- 
tion presents  the  conclusions  of  the  computations  discussed  here 
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2.   THE  MODEL 

We  consider  a  straight  cylindrical  plasma  column  whose 
equilibrium  profiles  are  considered  to  have  variation  only  in 
the  radial  direction.   The  plasma  is  considered  incompressible 
and  has  a  density  profile,  p  (r) ,  that  extends  to  a  plasma- 
vacuum  interface  where  the  density  is  zero.   The  equilibrium 
magnetic  field,  B  (r)  has  both  longitudinal  and  azimuthal  com- 
ponents.  Finally,  the  column  is  subjected  to  a  rotation  with 
rotational  frequency  Q(r). 

2.1   Basic  equations 

To  study  this  system  we  employ  the  following  mathemical 
model : 

||-  +  V-  (pv)  =  0  (1) 

3y 
P  ^rr  +  PY'Vv  =  -  Vp  +  JxB  (2) 


8B 

j^   +  VxE  =  0  (3) 


J  =  —  VxB  (4) 

-   y0   - 


E  =  -  V*B  +  nJ  +  -  JXB  (5) 

~  ~     -    p  -  - 


V«v  =  0  (6) 


V-B  -  0  (7) 
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where  v  is  the  plasma  velocity,  B  the  magnetic 
field,  p  the  plasma  pressure,  J  the  current  density, 

E  the  electric  field  and  n,  e  are  resistivity  and  Hall 
current  parameters  respectively  which  are  taken  to  be 
constant.   The  mathematical  model  specified  by  Equations 
(1)  -  (7)  differs  from  Ideal  MHD  only  through  Ohm's 
Law  which  includes  the  effect  of  finite  resistivity  and 
the  Hall  effect.   Our  aim  in  this  report  is  to  study 
the  time  behavior  of  small  perturbations  about  an 
equilibrium  specified  by  v  (r) ,  B  (r) ,  p  (r)  and  p  (r) . 
The  perturbed  velocity  field,  magnetic  field,  density, 
pressure,  current  density  and  electric  field  are  repre- 
sented by  u,  b,  6p,  6p,  j  and  6E  respectively. 

The  linearized  version  of  equations  of  (1)  -  (7) 
are: 


|t  fiP  +  V*(<5PY0)  +  v'(P0y)  =  0  (la) 

tt  p   u  =  -7p*  +  jxb   +  J  xb  -  p   u*Vv     (2a) 

dt     kO  ~       ^      ~  -O     -O  -      O  ~    ~0     v    ' 

-  PQ  Y0'Vy  -5pyo.Vyo  +  A_7(B-b) 
It  b  =  -V*<5E  (3a) 

j  =  i-  Vxb  (4a) 

o 

6E  =  -uxb   -  v  xb  +  nj  +  —  {jxB   +  J  xb  -  ^  J   B  }  (5a) 
~  ~o    ~o  -     i         p    i    -o         ~o  -    p   -o  ~o 

Ko  Ko 
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V-u  =  0  (6a) 

V-b  =  0  (7a) 

The  starting  set  of  equations  for  the  computation 
can  be  obtained  from  system  (la)  -  (7a)  by  Fourier 
transforming  the  perturbed  quantities  in  the  ignorable 
0,  z  directions  and  time,  t.   As  mentioned,  we  consider 
only  the  mode  m  =  0,  k  ^  0,    where  m  is  the  azimuthal 
wave  number  and  k  is  the  axial  wave  number.   Thus, 


b  =  b(r)  exp(iLu3t  -  kzj  ) 


u  =  u(r)  exp(i[ut  -  kz]  ) 


etc. 


For  ease  of  notation  we  will  omit  superscripts  from 
the  perturbed  transformed  variables. 

After  transforming  equations  (la)  -  (7a)  and  rear- 
ranging we  arrive  at  the  following  sixth  order  system  of 
ordinary  differential  equations: 

<*  f„   2  1  diT^\  •  f   d   ,1  d(rbr\    v2o  fr02  dPo 

dF(pou   F  dr-}  "  x  V~   dF  (F  -^r— )=  ~k  po(ril   ^dF 

-J    +    4Q2    +    2rfi    ^)5      +    i(i    [-    d!|  +    1   df    .    k2f2 

dr    ^r  p  ,    2         r   dr  J 

o        dr 

_    2k^  2k2f    .         2B9    k2         . 

w       Jo,z)br  -    <~&r  n  +  T"        )      b6  (8) 
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n   d   /Id   .  ,  .   _    r.  efk   , 

IT  dF  (F  dF  (rbr})  -  "  £^r   +   W   b6 
o  Ko  o 


,2 

+  (iw  +  £_E1  +  —  k  J    )  b  (9) 

Uo    pQ    o,z'   r 


£f   d   .Id  nk  d   .1  d   , 

Tr  (-  xr  (rb  ))  -  —  j-  (-3-  (rbQ)) 


Vi  P_  dr   r  dr     r      u   dr   r  dr 
o  o  o 

nk2    if2       irk   2Bfl    Bfl   dpo 

p„    up  p      u    rp       2   dr     6 
o     o  o     Ko    o    p 

o 

-  ku[J   „ ^-]  +  2kfft  +  krf  ^)  C 

o,z     y  r  dr    r 


,  ,e   .   d   ,  1  df.    k2f\    .   dfl   kf   ,   ,  . 

+  (IT  Cr  dF  (F70  dF>  +  "IT1  +  kr  dF  ~  ^o  Jo,z)  br    ao, 


where  f  =  kB    .   Furthermore,  we  have  considered  the  equilib- 
o,  z 

rium  velocity  field  to  be  a  rotation,  v   =  rft(r)6  where  9 

~o 

is  a  unit  vector  in  the  azimuthal  direction. 
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2.2   The  equilibrium 

We  consider  an  equilibrium  magnetic  field  profile 

with  vanishing  B     on  a  cylindrical  surface,  r  =  r  ,  in  the 

o,  z       J  o 

plasma,  i.e.,  a  reversed  field,  while  the  component,  B  , 


o'e 


of  the  field  is  finite  on  r  =  rQ. 


Since  the  plasma  is  incompressible,  we 
may  impose  an  arbitrary  density  profile.   Finally,  as  men- 
tioned, the  equilibrium  flow  is  taken  to  be  a  rotation  with 
variation  in  the  radial  direction. 

The  above  equilibrium  may  be  summarized  as: 


B  (r)  =  B^  Q  9  +  B     z  ,    B    (r  ) 
-o       o,9       o, z         o, z   o 


v  (r)  =  rft  (r)  6 


p  (r)    arbitrary 


2.3   Dimensionalization 

We  dimensionalize  with  respect  to  the  equilibrium 

values  at  the  neutral  surface  r  =  r   as  follows: 

o 


n      -  _   e    df  .   , 
o —   /   £  -  -5—  (r  ) 

2  d  Ly    u      dr    o 

cj  a  y  Ho  oMo 

o    o 
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o     -  P°(r)    f  -        *  df  _   dF 

Po  ■  PQ(ro)   '      adflro)   '  dr  "  df 

dr  dr   o 


(*   ) 


,*2^             aJI       °  r-r 

d  f  ,   .      dr  o    —   r 

— =■  (r  )  =  a  — - t-e =  a,  ,   t  =  — —  ,  r  =  — 

dr2    °       *l(r  )  X         a         a 
dr   o 


dp  B  n 

Ds  "    po  ^F   '   B8  (ro)  "  72~7" 

o  a   dr  .   , 

ar(ro) 


z     df  ,   .  '       u  ] 

adF(ro) 


tt   _  2fik    t-        ,  dQ.        2ak  dft 

u  dr     w   dr 

o  o 


a      df  ,   .  lbr 

o  o  dr   o 


qfl  = 


ikb0 


8     2  df,   , 
dr   o 


With  these  dimensionalizations 


f (t)  =  t  +  -1  t^  + 


and, 


J     =  a_x  +  ... 
o,z     2 
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We  rewrite  equations  (8)  -  (10)  in  dimensionless 

form  for  a  layer  centered  about  r   that  is  sufficiently 

o 

small  to  neglect  variations  in  the  equilibrium  quantities. 
This  layer  will  be  referred  to  as  the  boundary  layer, 
and  the  equations  governing  the  behavior  of  the  resistive 
layer  will  be  known  as  the  resistive  boundary  layer  equations 
or  RLE.   To  leading  order  in  x  and  omitting  the  super- 
scripts on  the  dimensionless  quantities,  the  RLE  are: 


2  2 

-d£  dq  -,     rD  r       ,n             -,-, 

2     r         ^r  ,n2ri  o  s,  o0dP.     .  2  2.  c 

dx^  dr  4  ^  a 


i        o         if2ka  O      2iBfl 

-(-   -  k2T-a. 2-)q   +  (iil 2.)q.         (8a) 

r         1     y    r    y      r   ^9 


Y    r     y 
o  '        '       o 


,2 

d  q 


dx 


2^  =  -YTCr  +  ETqQ  +  (Y+k  H  +  ieka2T)q] 


(9a) 


d2qe          d2qr  ,      2  2 

dx2  -«    ^2  =     (Y    +    nk2    +    I-    +    iekBft[?    -    Ds])q( 


r^  ho  2lkBo 

+    (    Ds    [1+k    a2Bel    +   ^    [1   +    2§fl    +   Y[a2T    "^]Kr 

1   )c  Pi  T" 

"(eta,    -   ±-  +   Ds   +   k2x]    -  ~f    x2    +   ^  g)qr  dOa ) 

o 
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3.   ANALYSIS  OF  THE  RLE 

3.1   Scaled  equations 

We  wish  to  study  the  interesting  behavior  governed 
by  the  RLE  in  the  limit  of  small  resistivity  and  Hall 
parameter  n  and  e.   Following  earlier  work  (COPPI,  et  al., 
1965)  we  adopt  the  following  scaling  based  on  a  resistive 

layer  thickness  of  the  order  of  the  skin  depth  with  respect 

2 

to  the  diffusion  time  scale,  i.e.,  yr      ~  n  where  t  ~  y. 


Y  =  4>s/T=<J)X,r|  =  <+>n   ,n  =  0(1) 


£  =  <{)£ 


5r  =  £  +  ...»  qr  =  <H  + 


qQ  =  C  +  • • • 


Also  denote: 


2J1  dr   '        r 

o 


where  the  superscripted  variables  are  considered  to  be  0(1) 
quantities  while  4>  is  a  small  scale  factor.   The  scaling 
of  the  Hall  parameter  and  angular  velocity  were  chosen  so 
that  terms  with  e  and  Q   would  neither  dominate  nor  disappear 
from  the  scaled  equations  in  the  lowest  order. 
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Finally,  we  again  omit  the  superscripts  in  writing  the 
scaled  equations: 


s2  d!|  +  x  df|  =  Q2(1  .  fo  D   +  p)€  +  («|  .  ifi)c  (11) 

dx      dx 

,2 

n  =-l  =  -sx£  +  exc  +  s^  (12) 

dx 


^2.      ^2,  2  r  Dc 

n   dj;   _  £d_|  =  {s  +  x   +  ie6{1  _  Jpm 

dx       dx 


+  (x«[l  +  P]  -  i<5S)£  -  Pity  (13) 

Making  use  of  (12) ,  equations  (11)  and  (13)  can  be 
rewritten: 

S2^=  (*?S  +  fl2(1  _  !o^S  +  p])  c+  (Q|  _  ex! 
dx2     n  4  b     n 

-  15)  c  -  —  i>  (Ha) 

2  2      2  r 

n  -^   =  (S[l  +  |-]  +  ^-  x2  +  ie5[l  -  -|  Ds])c 
dx 

2 
+  (xfi[l+P]  -  i6S  -  ^*-  S)5  +  (—  -  P")^  •  (l3a) 

Equations  (11)  -  (13)  show  that  in  the  presence  of 
azimuthal  components,  Bfl ,  of  the  equilibrium  magnetic  field, 
overstable  modes  are  possible  due  to  the  imaginary  terms 
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in  the  equations.   For  this  analysis  we  suppress  these 
terms  by  considering  only  long  wavelength  modes  with 
sufficiently  small  k  which  results  in  the  elimination 
of  6  and  a   from  equations  (11)  -  (13)  to  leading 
order. 
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3.2   Asymptotic  behavior  of  RLE 

The  system  (11a) ,  (12) ,  (13a)  constitutes  a  sixth 
order  homogeneous  linear  differential  equation  with  an 
irregular  singularity  at  x  =  °°.   The  system  may  be 
rewritten  in  a  standard  form  for  irregular  perturbation 
analysis  (WASOW,  1965). 


du 
dx 


=  x  A(x)  u 


which  has  a  set  of  solutions  of  the  form: 


A  x         D 
u  =  c  e  x   z  (x) 


where  c  is  a  constant  matrix  and  A  ,  D  are  diagonal 
matrices . 

The  elements  of  A  determine  the  exponential  nature 
of  the  asymptotic  solution  to  the  inner  equations.   These 
elements,  X.,  are  the  eigenvalues  of  the  matrix, 


A(»)  = 


0 
0 
0 

1_ 

nQ 


o 
o 
o 

E 


nQ' 


—  (h   +  — ) 
n  nQ    n 


o 
o 
o 


o 

0 


1 

0 
0 


0 
0 


0 

1 

0 


0 

0 


0 
0 

1 


0 
0 
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Thus, 


Al,2  =  °'° 


2  1/2 


A3  4  =  -±72     (l  ±     (1  +  4~}    } 

2  !/2 

Xs  6  = Tn  (l  *  (1  +  f-)      } 

b'b     (ns)1/z  l  4 

u        2     £2S 

where  r   =  . 

n 

At  this  point  it  is  not  clear  whether  A  , 
k  =  3,4,5,6  have  pure  or  oscillating  spatial  exponential 
growth  or  decay,  or  whether  they  have  pure  oscillatory 
behavior  as  x  -*■   °°.   In  fact,  we  can  shovv  that  only  if, 


S  =  3 (-1  +  -) 

/  2.L.1  x  2      "  a 
(a  +l)e 


for  a  an  arbitrary  real  number,  can  A   be  pure 
imaginary  in  which  case 


A,  =  ±  -  i  . 
k      n 

Thus,  all  these  solutions  either  decay  or  grow 
exponentially  as  x  ->■  °°  in  the  class  of  solutions  with 
positive  growth  rates,  i.e.,  Re  S  >  0  , 

Finally  we  consider  the  two  algebraic  equations 
corresponding  to  A   and  A  .   To  determine  the  asymptotic 
behavior  of  these  solutions  as  x  -»■  °°,  we  consider  the 
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dominant  terms  from  equations  (Ha)  and  (13a),  as 


sc  e  n  -   sip 
n    n    xn 


(a-s£)5+(i  +  £i)c:^    , 


Solving  to  leading  order  in  x  for  £  and  y  yields 


i  :  i  ,  d4) 


■^§   *    •  (15) 

x 


Thus, 


■f?       •  (16) 


Also    from    (12) , 


dj^   =    d2(xg)  d2       ,3gr  (17) 

^2  ,2  e,2ls"J 

dx  dx  dx 


Replacing  (14)  and  (15)  in  (11)  and  noting  that 


'  £  <*>  ■  fe  <*2  i' 
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yields,  to  leading  order  as  x  -*■   °°  , 


§_  (x2|i,  .  .Q2  £D   -m  ,18) 


which  has  the  solutions, 


€  -  xh  ,  -"  "h 


x 


r  D 
where  h  =  -  j  +  i(l  -  ft2  (-j-^-  -  P))1/2.   In  order  to  arrive 

at  (16)  we  made  use  of  the  fact  that  in  view  of  (14), 

h-1     -2-h 
Y  -  x     ,  x 

In  summary,  the  asymptotic  solutions  to  system 

(11)-  (13)  take  the  form,  as  x  ->  <*>    , 

2    ^    2 


-1     r         r2  1/2   x2 
exp  { *    (|  ±  (1  +  |_)    )  |-} 

(ns)  7   ^       4      ^ 


h     -h-l 

X    ,  X 


2  r 

where  r  =  £-§•  ,  h  =  -^  +  i  (1  -  fi2  (-^  D   -  P)  )  1/2 
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4.   SOLUTIONS  OUTSIDE  THE  RESISTIVE  LAYER 

The  behavior  of  the  plasma  outside  a  narrow  region 
surrounding  the  neutral  surface  is  governed  to  leading 
order  by  ideal  MHD,  i.e.,  e  and  n  can  be  taken  to  be  zero. 
In  terms  of  the  same  dimensionalizations  as  the  resistive 
layer  equations,  the  outer  equations  reduce,  near  marginal 
stability,  i.e.  Y  =  0 ,  to  the  following  second  order  differ- 
ential equation: 


d    ,f2    d     ,    _    .,    f2,2r  Po      .2,_      2   d«,  _  ,       . 

dF(F"  dF(rV)-f   k    ^r  =   "  -T  rfi    (Dsl  d7Kr  (l9) 


Near   r      equations    (11)  -     (13)    reduce    asymptotically    to 


d       ,    2   dCr.  n2  ,roDS 

dT    {T      dT]    =   "fl    ("T-  "   P)?r  (20) 


which  has  solutions, 


t-  h      -1-h  ,_ ,  , 

£  =  T     /    x  (21) 


where  h  =  -|+|(l-  fl2  (rQDS  -  P))1/2.  (21a) 
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For  the  purposes  of  numerical  integration,  it  is  con- 

venient  to  set   W  =  -r= —   which  transforms  (19)  to, 

as 
r 

dx 
|£  +  P(x)  W2  +  q(x)  W  =  1  (22) 


W(0)  =  W(l)  =  0 


where, 


-2 
,  .    1     2  df/dx    .2    DsX" 

P(x)  =  "2  "   xf +  k  "  —7T- 

x  4f 


and, 


,  .      .2  df/dx    1. 
q(x)  =  -  (   — | +  -) 


df  .. 


From  (22)  and  the  asymptotic  expressions  f  ~t  ,  ^   ~    1 , 
we  find  that  the  criteria  for  entering  the  asymptotic 
regime  governed  by  (20)  is  that 

M  -  |   f1  _2  I  «  1  - 

In  this  formula,  the  value  of  x  such  that  M=.  01  may  be 
taken  as  a  measure  of  the  thickness  of  the  boundary  layer,  i.e., 
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Tb.£.    '01   -T-  (23) 


We  note  that  according  to  (23)  the  boundary  layer 
thickens  as  the  plasma  rotates  faster. 

5.   MATCHING  INNER  AND  OUTER  SOLUTIONS 
In  order  to  match  the  outer  solutions,  the  two 

exponentially  growing  solutions  to  the  inner  equations  must 

be  set  equal  to  zero,  leaving  the  two  algebraic  solutions 

for  matching. 

The  condition  that  the  algebraic  inner  solutions  match 

the  outer  ideal  MHD  solutions  can  be  expressed  as: 


outer       inner 

;r(T)  -  n 

T+0  X-*"=° 


lira  £r  (T)  -  lim  £  (x)  (24) 


From  (18)  and  (19)  it  is  clear  that  the  inner  solutions 
can  be  matched  asymptotically  to  the  outer  solutions  in  the 
sense  of  (24) . 

We  will  be  investigating  a  class  of  resistive  eigenmodes 
called  tearing  modes,  that  are  sensitive  to  the  matching 
conditions  between  the  inner  and  outer  regions.   We  subdivide 
the  plasma  radius  into  three  regions.   Region  I  is  the  outer 
region  surrounding  the  magnetic  axis.   Generally,  when  treat- 
ing boundary  layers,  the  boundary  layer  is  called  the  inner 
region,  and  the  rest  is  called  the  outer  region.   Here,  how- 
ever, we  have  two  outer  regions,  one  of  which  is  in  some 
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sense  inside  the  boundary  layer.   We  can  just  call  it  region 
I.   Region  II  is  then  the  inner,  resistive  boundary  layer, 
and  region  III  the  other  outer  region. 

From  Equation  (21)  we  have  seen  that  the  solution  in 
Region  I  near  the  boundary  layer  approaches  the  form: 


^r,I^AI  I^+Bjlxr1-* 


(25) 


Likewise  in  Region  III, 


«r.Hl  =  *IIIlTlh  +  BlllM"1"h 


(26) 


Similarly  in  Region  II, 


C  z 


Aii,LWh  +  "illM"1"*  -  *  *  - 


AII.RlX|h+BII,RlXl"1"haSX* 


The  matching  condition,  (24)  can  be  expressed  as, 


B. 


,1.  2h+1  ^L 
V        BII,L 


(27) 


and, 


III 


B 


III 


1   2h+1  AII,R 
*        BII,R 


(28) 
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for  2h+l  =  {1  -  Q2(r  D„  -  P))1//2  >  0  and  where  <J)  is  the 

o  s 

small  scale  factor. 

A       A 
We  expect  for  most  cases  that  _I  and   III  will  be 

BI      BIII 
finite  numbers.   Thus  to  obtain  matching,  the  growth 
rate  must  be  adjusted  in  the  resistive  layer  so  that 
ATT    -  0.   The  exceptional  case  occurs  when  B   vanishes. 

X  X  /  Li  I 

According  to  Newcomb's  stability  criterion  (NEWCOKB,  1960)  this 
is  the  condition  for  marginal  ideal  stability.   It  is  not 
surprising  that  a  system  that  is  close  to  being  unstable 
in  the  ideal  theory  will  have  a  strongly  modified  resistive 
instability. 

These  matching  constants  can  be  calculated  numerically 
from  the  outer  layer  equations  as  the  limit 


d   ,  l+hr. 
AI      AIII    lim   ,-  -4h-2  dx  ^X    K\  ,   . 

it    -    bttt  =  x-+-  (+x — )        (29) 


d_   h 
=  lim   .-  4h-2  dx  ^x  v; 

l+X 


x-»-+cd         d   ,  -1-h,. 
dx"  (x    *> 
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(30a) 


In   terms   of   W,    (29) can  be   rewritten, 


Ani    iim  -  iTr2h_1    h+i  +  5 
BIII    T-°  I  — — -  1 

_h  +  W 


A!    i •     I  |-2h-l    ...    T 
=-  =  lim  -  | t |         h+1  "  w 
I    t-0  [ ] 

-h  -  s 

w 


The  quantity  matched  to  the  inner  layer  is, 


A'  m   J4li  _  ^i  m  (3(JbJ 

BIII    BI 


Our  objective  in  this  report  is  to  gain  some  infor- 
mation about  the  eigenvalue  problem  by  purely  analytic 
means.   To  do  this,  we  consider   the  growth  rate, 
S,  and  h  -  0.   It  will  prove  convenient  to  match  the 
r-component  of  the  perturbed  magnetic  field  at  the  edge 
of  the  resistive  layer  rather  than  the  plasma  displacements, 
in  which  case  the  asymptotic  solutions  are 


q     -  B   +  Ax   ,   as  x  ->•  0 

1     /    L  X  X 


*r,IU  =  Bin  +  AinT   '   aS  T  *  ° 
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and  , 


i>  - 


BII,L  +  AII,L  X 


as  x-»— c 


BII,R  +  AII,R  X     as  x  +  ~ 


As  in  (27)  and  (28)  the  matching  conditions  become, 


I  AlI'L      ,      itlls*  =    o(l) 


<t>    B 


II, L 


<p    B 


(30) 


II, R 


We  have  seen  that  the  behavior  of  the  exponential 
solutions  of  the  RLE  are 


£  ~  exp  { ± 


kx' 


2(nS) 


172 


}   for  k  constant 


as  x  -*±».   Thus  as  S  becomes  small,  we  must  consider  the 

1/4 
length  scale  to  become  small  as  S    ,  otherwise  the 

exponential  solutions  would  grow  infinitely  rapidly  making 

the  equations  badly  posed.   That  is,  there  is  an  inner 

layer  of  significant  behavior  that  is  even  thinner  than  the 

1/4 

resistive  skin  depth   by  another  factor  of  S 

We  find  the  interesting  behavior  within  this  new  scaling 
by  expanding  the  parameters  and  variables  of  equations 
(11)  -  (13)  in  powers  of  S  as  follows: 
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5  ■  5o  + 


*  =  s1'4  *o  +  s7'4  ♦  ♦ 


-  s3/2    + 


^O 


x  =  sl/4  -  ,    e  =  ~j^   e  ,  n  =  s3/4  n 


Equations  (11)  -  (13)  reduce  in  the  zero-th  order  to: 


a2i^ 
n  — =2  =  0  (31a) 

dx 

2 


d2?o      2      r2  -  "2 

— «2  =  xz    [i+-  ]z,     +    (xflti+p]    -  ^-)£     +    (e  -  -  PW 
dx2  no  no  n  o 


In  the  first  order: 


Thus 


n  ~   =  -x£   +  ex^   +  *  (31) 

,  -2        O       o     O 
d  x 


(0)     (0) 


*o  =  BII   +  AII  X 
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But  since  A 
(0) 


II 


0  for  matching,  we  are  justified  in 


setting  A     =  0  in  the  zeroth  order  and  accounting  for 
the  x  dependence  in  the  next  order.   In  fact,  from  (31), 
and  in  view  of  (14)  and  (15),  we  have 


d2ii 


dx 


1  =  0 


as 


x  -»•  ±°° 


which  has  solutions, 


*1  = 


(1) 
A  T    x  as  x->- c 


(1) 
AII,R  *  as  x  -> 


where  the  constant  solution  has  been  incorporated  in  the 
zeroth  order  solution.   Thus  the  total  ty    solution  is, 


7/4   (I)-     1/4  (0) 


And  since, 


If     (x>  "  "IT?  ^  <*' 

aX  Sl/       dx 


making  use  of  (25),  (26)  and  (31),  we  arrive  at  the  dis- 
persion relation, 


!a    a' 


(32) 


where, 
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A   = 


B 


III 
III 


BT 


(33) 


and, 


** 


lim  dxU)  dx(X  ' 


S5/4  lim 


x"~  <t>    6  S1/4    3  ♦  X+-  J 


-X 


-X 


[-X^o  +  £  X  SQ 


+  B]dx 


(34) 


where  we  have  set  0  =  B„'E  tb  . 

II     o 

Thus  the  condition  that  a'  =  0(1)  requires  that 


*  =  s5/4 


(35) 
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6.   ANALYSIS  OF  TEARING  MODE  EQUATIONS 

6.1   The  Scaled  equations 

From  the  original  scalings  of  section  5,  and  equation 
(34),   the  tearing  mode  scalings  can  be  rewritten  in  terms 
of  the  scale  factor  9  as 


„  .  ft  ,   Y  .  e3o  ,  §  -   e§A 


t  =  62x  ,  e  =  9e  ,  a   =    9A 


qr  =  e2^o  +  e4^1  +  . 


C  =  5  +  .. 

r 


(36) 


^9  =  6  ^o  + 


*Q  =  B 


and   letting  P  =  ^?  ^ 

2A  dr 


£D, 


and  q 


(36a) 


The  following  tearing  mode  equations  result  from  scaling 


the  basic  equations  (11)  -  (13)   according  to  (36), 


omitting  superscripts  and  considering  |p  - 


o  S 


<<  1 


<J25. 


dx' 


.2 

(\  (1 
Q 


o  S 


+  P]  + 


2 

Qn  ^o 


+  ( 


Ax 


ex 


Q  n 


>£. 


1_ 

r 


+  (- 


-  — )  B 
nQ 


(37) 
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d2C 


dx 


o.    (2i[1  +  p+ql    -4oKo+    4  +  ^<o 


,   e  r    1  t   A ,         .         exQ 

+  (~  n[ai"F-]  "nCp+q]  +  "2" 
o  n 


)B 


d2ij; 


dx 


I-  -  2fi  c  +  £21  ,  +  Q  B 
2     n   o   n   o   n 


(38) 


(39) 


where  the  superscripts  have  been  omitted  from  the  scaled  variables, 
6.2   Asymptotics 

We  study  the  asymptotics  of  system  (37)  -  (39)  by 
expanding  the  independent  variables  as  follows: 


2 
x 


5  =  =±  +  -4 


?1    ?2 

X 


(40) 


(41) 


Replacing  these  expansions  in  (37)  -  (39),  we  obtain 
the  result  that 


?1  =  °   '  ?1  =  B  ' 


C2  =  -8QA  ,  S2  =  [g  (~  ~   a1)  -  cA]B 

o 


(42) 


and, 

.  r  D 

t|'1  =  6[1  +  1-5-  In  x  +..-]  +  [c   +  B(a  -  ~)  In  x]x  +, 

o 


where  B  and  c   are  chosen  to  match  the  outer  solution 
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6.3    The  dispersion  relation 

In  order  to  study  system  (37)  -  (39)  analytically 
we  are  motivated  by  equations  (42)  to  make  the  following 
change  of  variables: 


z  = 


X 


(Qn) 


174 


5  = 


B 


(Qn  ) 


1/4 


We  also  define  the  following  parameters, 


2    e  Q 
r  =  ^—^ 


and  q  =  P  +  q 


A  = 


,3/4 
1/4 


(43) 


With  these  substitutions  (37)  -  (39)  can  be 

r  D 
rewritten,  assuming  that  |  - u.         -    P|  <<  1: 


d27     -?      r  Dq      ?  9 

— I  =  [A  (1 J-^)  +  7.    \\    +  A[Az  -  rz^H  +   (a1 

dz 


l     V4 

i-)  H 2 

rJ   5/4 

O  Q 


(44) 
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dz2 


-  (zA(l  +  q)  -  rz2K  +  z2d  +  r2K 
A 


+  ( 


1/4 
rz    en  '   /     1  >    - 

r '  ^74^  (ai  -  ? )  -  q 

A     Q    A         o 


(45) 


d2!    Q^ 

dz2   n1/2 


(1  -  z£  +  rA?) 


(46) 


In  the  text  that  follows  we  will  simplify  the 

1 


computation  by  taking  a,  = 


The  dispersion  relation  (32)  -  (34)  can  also  be  rewritten 
in  terms  of  these  substitutions: 


where, 


ifeiT-  &L  F 

Sl_dxj   n3/4 


f  z 


F   =    lim 

2-»-oo     J 


(1    -    z£    +    rA?)dz 


-z 


(47) 


and, 


lfd^l      1    1 
3LdxJ      3   X 


im    r^i    (x)_    d± 

+°°    ldx    lx;       dx    (    X}  J 


=    A, 


Thus,  the  dispersion  relation  is 


A'  =  A 


(48) 
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and  making  use  of  (36)  we  have: 

-  .  -3/5  (A',  4/5  ,49) 

where  the  dimensionless  quantities  y  and  n  have  been  used 
in  place  of  the  scaled  variables  Q  and  n  /  and  F  depends 
on  the  parameters  in  equations  (44)  -  (46)/  namely, 
q,r,A,  and  D  • 

6.4   A  special  case 

Our  objective  here  is  to  find  a  special 
class  of  analytic  solutions  to  the  inner  layer  equations 
(44)  and  (45)/  and  use  these  solutions  to  evaluate  the 
dispersion  relation  (48)  explicitly. 

The  case  we  consider  is  characterized  by, 

r  D 

~-^-  «    1,   r  <<  1,   P  -  0,   q  =  -  1,   A  -  0(1) 

This  corresponds  to  a  special  case  of  a  theta  pinch 
equilibrium  with  rigid  body  rotation  and  density  gradient 
opposite  in  sign  to  the  angular  velocity. 
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From  (44)  we  see  that  only  the  odd  part  of  E, 
and  c,    contribute  to  the  dispersion  relation.   Furthermore, 
in  the  limit  of  small  Q,  the  second  term  in  (44)  can 
be  treated  as  a  perturbation,  thus  only  the  odd  part 
of  %   must  be  determined  in  order  to  derive  the  dispersion 
relation.   But,  referring  to  (43)  and  (44),  and  using 
the  result  of  the  asymptotics  summarized  in  (14)  and  (15), 
all  terms  involving  r  may  be  treated  as  pertubations  and 
neglected  to  leading  order.   Thus  (44)  -  (46)  reduces  to: 


d2C     2-    -2    -  72    - 

— f  -z^-AF   =    zh      t,   ~    z 

dz 


(50) 


&...*,. 


dz- 


z-  c  =  1 


(51) 


6  ldxJ    ,  _1 


(nQ) 


/4 


¥1 

1/2 


(1  -  zC)  dz 


(52) 


where  we  denote  the  value  of  A  corresponding  to  the 
special  case  q  =  -1  by  A  .   We  let, 


I   -    e"z  /2  I    b2k  H2k(z) 


(53) 


l  «  rz  /2 


^  a2k+l  H2k+!(Z) 


(54) 


wh 


ere  H_,  ,  H_,  ,n  are  Hermite  polynomials. 
2k     2k+l  r 
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and  make  use  of  the  identities, 


1  =  /2 


=  /2  e 


-z2/2 

H   (z) 

2k 

-z2/2 

9  V 

2   r  (k+l) 
H2k+l(z) 

22kr (k+l) 

Solving  for  £  we  find, 


[2  +  A  ( 


o v4k+l 


4k+5 


)] 


2k+l 


22k+2r(i+i) 


(4k+3  +  A2 
o 


Replacing  E,    in  (51)  ,  and  using  the  orthogonality 
relations  for  Hermite  polynomials  yields, 


-5/4 
A  =  ^74-F(Ao) 

2n 


where, 

-2 
3+A^  . 

-2   r(-l^     -2  T{V 
F(Ao)  =  -4(1  +  A2)  -— ig-  +  2A2  -f- 


5+A' 


nf) 


Growth  rates  can  be  obtained  by  solving  the  dispersion 
relation,  A   =  A,  for  the  growth  rate  as  given  by  (49) 
where  the  quantity  A*  is  determined  by  integrating  the  ideal 
equations  of  motion  forward  from  r  =  0  and  backward  from 
r  =  rw(rw  =  wall  position)  and  then  forming  the  logarithmic 
derivative  as  indicated  in  Eq .  (30a).  From  (49)  we  notice 
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that  the  stability  criterion  is  A '  >  0  ,  that  is,  for  positive 
values  of  A   the  plasma  is  unstable  and  positive  growth  rates 
can  be  calculated.   For  A   <_  0  ,  which  depends  only  on  equilib- 
brium  quantities  outside  the  resistive  layer,  the  plasma  is 
stable.   We  note  that,  in  general,  the  dispersion  relation 
will  have  to  be  solved  iteratively  since  from  (43),  A   depends 
on  the  growth  rate. 

6.5   Validation  of  the  analysis 

The  ultimate  goal  in  this  work  is  to  compute  the  stability 
(growth  rates  and  stability  boundaries)  of  realistic  plasma 
configurations  to  the  tearing  modes  considered  here.   Before 
proceeding,  it  is  improtant  to  validate  the  analysis  with  a 
specific  example.   Let  us  consider  the  special  case  of  the 
last   section.   We  will  show  by  a  numerical  integration  of 
the  governing  equations  in  the  resistive  layer  that  i)  the 
asymptotics  are  correct  and  ii)  the  analytic  solutions  given 
by  (53)  and  (54)  are  correct.   We  use  a  technique  (FISCHER, 
1969)  for  solving  Eqs.  (50)  -  (52)  which  converts  the  bound- 
ary value  problem  into  an  initial  value  problem.   This  tech- 
nique has  been  used  successfully  elsewhere  (JARDIN)  in  a 

number  of  resistive  boundary  layer  calculations.   We  have 

-2 

chosen  the  special  case  q  =  -1,A   =  0.75  and  have  numerically 

integrated  the  resistive  layer  equations  for  \    and  C  as  a 
function  of  z  for  0  <  z  <  4.   Shown  in  figure  la  is  a  curve 
of  \    vs.  z ;  the  crosses  denote  the  solution  as  given  by  the 
analytic  expression  (54)   and  the  circles  denote  the  results 
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from  the  numerical  integration  of  Equation  (50)  and  (51) 
Figure  lb  presents  a  similar  comparison  between  analytic 
and  numerical  solutions  for  If  vs.  z.   The  asymptotic 

expansions  for  c,    and  £  for  large  z,  i.e.  f  -.  1/z  and 

—       2 

£  ~  1/z  ,  are  borne  out  by  both  the  numerical  and 

analytical  solutions.   Further,  comparison  of  the  solu- 
tions for  the  entire  z  range  indicates  excellent  agree- 
ment.  We  believe  therefore  that  the  numerical  technique 
for  solving  the  RLE  equations  is  satisfactory. 
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7.   GENERAL  RESULTS 

7.1   Effect  of  Hall  current  on  the  growth  rate 

From  equation  (49)  we  see  that  the  effect  of  Hall 
current  on  the  growth  rate  is  entirely  contained  within 
the  quantity  F(A,q,r,D  ).   The  value  of  F  is  governed, 
for  the  most  part,  by  the  parameters  of  the  inner  layer, 
namely,  e,n  and  the  local  values  of  the  angular  velocity, 
ft,  and  the  density  gradient,  D  .   Thus  the  quantity  F 
may  be  considered  as  a  criteria  of  merit  which  gives  us 
a  means  of  evaluating  the  effect  of  Hall  current  on  the 
growth  rate  of  the  tearing  mode  as  compared  to  growth 
rates  in  the  absence  of  Hall  current.   For  example  in 
the  absence  of  Hall  current  and  the  case  of  rigid  body 
rotation,  e  =  q  =  0  and  F  =  F(A,0,0,D  ).   Thus  larger 
values  of  F  in  the  presence  of  Hall  current  compared  to 
no  Hall  current  would  according  to  (49)  tend  to  depress 
the  value  of  the  growth  rate  whereas  smaller  values  of  F 
would  tend  to  have  the  reverse  effect  of  enhancing  the 
growth  of  the  tearing  mode. 

It  is  important,  therefore,  to  examine  the  criteria 
of  merit,  F,  more  closely.   We  will  demonstrate  later  that 

for  realistic  values  of  the  parameters,  F  depends  only  on 

_->      _  -2  - 

A   and  q.   In  other  words  F  =  F(A  ,q) ,  and  the  parameter  r 

along  with  the  explicit  dependence  of  F  on  Dg  through  the 

r  D 
term   ?  s  in  equation  (44)  may  be  considered  as  small  pertur- 
4 

bations.   This  simplifies  the  analysis  of  the  inner  layer  con- 
siderably and  enables  us  to  illustrate  the  dependence  of 
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-2 

F  on  A   and  q  by  Figure   2.   If  we  consider  rigid  body 

rotations,  the  curve  q  =  0  corresponds  to  the  case  of  no 
Hall  current.   From  this  figure  it  is  evident  that  values  of 
q  >  0  result  in  enhanced  stability  to  tearing  modes  whereas 

values  of  q  <  0  result  greater  instability.   For  the  case  of 

eD 
—      s 
rigid  body  rotation,  from  (  43)  ,  q  =  .   Thus  rotation  of 

n 

the  plasma  in  one  direction  relative  to  the  value  of  the 

density  gradient  at  the  neutral  point,  enhances  the  stability 

to  tearing  modes  while  rotation  in  the  opposite  direction 

results  in  greater  instability.   This  result  was  observed  by 

KALEK  (19  72). 

Figure  2   also  shows  that  the  effect  of  rotation  on  F, 

and  therefore  the  growth  rates,  is  comparatively  small  for 

the  case  of  no  Hall  current  due  to  the  relative  constancy 

of  the  curve  q  =  0  =  e  compared  to  curves  with  larger  values 

of  q. 

It  should       be  noted  that  in  the  case  of  no  Hall 

r        - 

current   but   a    sheared   velocity   profile,    q   =   — °    -^    , 

2Q   dr 

Figure  2  shows  that  sheared  velocity  can  also  effect  the 

stability  of  tearing  modes.   This  point  will  not  be  pursued 
in  this  paper.   Hereafter,  we  will  specialize  our  analysis 
to  rigid  body  rotations. 
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7.2   Effect  of  the  equilibrium  and  position  of  the  outer  wall 
on  the  tearing  mode 

The  equilibrium  determines  the  parameters  of  the  inner 
layer  and  thus  influences  the  tearing  modes  through  the 
criteria  of  merit,  F.   The  equilibrium  also  affects  the 
tearing  modes  through  the  outer  equations  via  the  parameter 
A  .   Furthermore,  since  by  equation  (49),  the  value  of  A 
is  extremely  sensitive  to  the  position  of  the  outer  wall 

surrounding  the  vacuum  region,  we  see  that  growth  rates  and 

-2 

thus  values  of  the  inner  parameter,  A   defined  by  (43),  are 

very  sensitive  to  the  outer  wall  position.   Thus  values  of 

-2 

A   determine  not  only  the  effect  of  rotation  but  also  of  the 

position  of  the  wall  on  the  value  of  the  criteria  of  merit, 
F. 

We  now  show  how  the  Hall  current,  rigid  body  rotation, 
and  wall  position  effect  the  growth  rates  of  tearing  modes 
for  a  family  of  rigid  rotor,  Theta  pinch  equilibria.   To  be 
exact,  we  will  determine  an  overestimate  of  the  ratio  of 
growth  rate  with  Hall  current  to  growth  rate  with  no  Hall 
current,  and  show  how  these  ratios  are  affected  by  angular 
velocity  for  fixed  positions  of  the  outer  wall.   The  next 
section  is  devoted  to  describing  the  equilibria. 
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8.   RESULTS  FOR  A  SPECIFIC  FAMILY  OF  EQUILIBRIA 

In  this  section  we  apply  the  theory  to  a  class  of  reversed 
field,  isothermal,  rigid-rotor  Theta  pinch  MHD  equilibria. 
We  assume  that  all  the  mass  and  rotational  energy  is  carried 
by  the  ions  and  that  the  family  of  equilibria  is  parametized 
by  the  angular  velocity,  fi,  and  characterized  by  a  constant 
number  of  particles,  N  .   The  equilibria  are  one  dimensional, 
and  although  their  domain  of  definition  is  from  the  plasma 
center  line,  r  =  0,  to  r  =  °°,  the  density  is  effectively  zero 
for  r  >_  a;  where  a  is  taken  to  be  the  plasma  radius.   Further- 
more all  equilibrium  quantities  are  constant  from  the  plasma 
edge  to  the  wall  of  a  vacuum  region,  r  ,  surrounding  the  plasma 

8.1   The  equilibrium  equations 

The  equilibrium  is  characterized  by  the  force  balance, 


aV  (p  +  w~ fi2)  =  "2rp  (55) 


and  the  equation  of  state, 


p  =  cp   ,   c  =  ^—  (56) 

i 


where    K    is    Boltzmann's    constant,    m.     is    the   mass    of    an   ion,    and 
T   is   the  plasma    ion   temperature. 
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The  number  of  particles  per  unit  plasma  column  length 
(ions),  N  ,  is  determined  by, 


NT  _  2tt 

Li     —  

m. 

1 


r3. 


p (r) rdr 


and  expressed  in  terms  of  a  characteristic  radius,  r0,  as 


N     tt    ,  .  2 

L  =  m~  P(o)r£ 

l 


(57) 


We  consider  the  following  ansatz, 


p (x)  =  p (o)  sech  ( 


x 


"  V 


(58) 


and, 


(5+ 


B   =  B   |S+(l-fi)  tanh( 

2      °° 


'   ^2 
x 


>V 


-  £ 


°>] 


(59) 


where,  all  lengths  are  normalized  to  the  plasma  radius, 


i.e.,  x  =  — 
a 


Replacing  (58)  and  (59)  in  (55)  and  (56)  and  character- 
izing the  ratio  of  the  thermal  to  the  magnetic  energy  by, 


P  (o)u 
8  =  2c  2.   ' 

B  2 


and  the  ratio  of  the  rotational  to  the  thermal  energy  by 


2  ? 
«.**  -   *g   B 
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yields, 


0*2  .  26  (l-o)  sech  E, 

X 


6  =  d-o)2  sech2C 


x2  m   x2  sech2C0 
p       1+tanh  £ 


The  equilibrium  is  therefore  entirely  specified  once 
the  physically  meaningful  variables,  B  ,  p  (o)  ,  T,  N   and  Q, 

00  Jj 

are  specified.   However,  for  our  purposes  it  is  more  con- 
venient to  specify,  B  ,  N    £  ,  6,  and  x0 . 

In  terms  of  the  latter  set  of  parameters,  the  plasma 
temperature,  rigid  body  rotion  and  the  neutral  point,  x  , 
of  the  plasma  are  given  by, 


,    2tt6  (1-6)  (1+tanh  £  ) 

JT  = —  (60) 

m.  y   NT 
l   o   L 


m.  c 

T  =  _L_ 
K 


where , 


c  =  n2a2x2z     (1-tanh  Cq)  (^-) 
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and, 


1/2      1 •  1/2 

xq  =  x£  (1-tanh  CQ)  '   (^0+^  ln(l-2Y))  /Z 


The  object  of  this  study  is  to  show  the  effect  of  the 
Hall  current  dominated  boundary  layer  on  the  growth  rates, 
y,  of  the  tearing  modes  expressed  by  the  dispersion  relation, 


-3/5 
Y  =  n 


F(A2,q) 


4/5 


(61) 


where,  once  again,  the  superscripted  variables  are  normalized 
with  respect  to  the  boundary  layer  quantities  and  the  variables 
q  and  A   governing  the  inner  equations  are  defined  by  (36a) , 
and  (43).   Values  of  the  mismatch,  A',  are  obtained  by 

integrating  the  outer  equations  as  described  in  section  7 

-  2  • 

while  F(A,q)  is  obtained  by  integrating  the  inner  equations 

as  described  in  section  8.   .   Since  y  occurs  on  both  sides 

of  the  equation, (61)  must  be  solved  iteratively. 

The  inner  variables  and  the  normalizations  for  (61)  are 

computed  as  follows: 


*  =  !?6  (Co  +  I   ln  (1"2Y)) 


(62) 


A2  " 


-3/2 
Y 


(63) 
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where, 


2 

-,        2e      6    cy 

2.  n< 

r     = 


nd-26) 


"  =  ifc^  (64) 


r   _    2fla 


kV 


k   =    2n  a 
L 

and,  (L  =  perturbation  wavelength) 

0         4a2TT[C      +   ^   ln(l-26)j  (1+tanh    £    )     (1-26) 
u      m.    NT 


and   from   equation    (21a) 


i  l  c*2       V2 


(65) 
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8.2   Results  of  tearing  mode  analysis 

We  carried  out  the  analysis  of  the  above  equilibrium 
for  the  following  range  of  parameters: 


Bm  =  1    '    ^  =  l    '    xo    =    1    ' 

<=  O  X 

.     no     ,.-8               MHO-KG                           ._    MHO-KG  ,,        .^-in-7 

e    =    1.02x10         e    -    — 2 ,    n    =    -12    —z ■,    u    =4ttx10 

m  weber  m 

2 

K  =  1.38xl0-22  kg  "  m2  ,   m.  =  3.2xl0~27kg 

o,  sec 
k 

a  =  .05  m  ,    —  =  .05 
NL  =   lxlO16  -  lxlO18  (ions) 


and, 


6  =  .025  -  .425 


where  6  is  the  normalized  angular  velocity  defined  by  (60). 

As  we  discussed  in  the  previous  sections,  the  effect 
of  the  resistive  layer  on  the  tearing  modes  is  contained 
entirely  in  the  criteria  of  merit,  F,  whereas  the  outer 
regions  affect  the  tearing  modes  primarily  through  the 
mismatch  A  .   Finally,  there  is  a  feedback  between  outer 
and  inner  regions  that  is  particularly  sensitive  to  the 
position  of  the  outer  wall,  since  the  position  of  the 

outer  wall  directly  affects  A   and  indirectly  affects  F 

-2 
through  A  .   As  we  mentioned  above,  F  depends  to  leading 

—  2 
order,  on  the  parameters  q  and  A  ,  since  for  the  equilibria 
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-3        r  D 
being  considered,  the  parameter  r  <  10   while   °  s  ,  accord- 
ing to  Table  1,  varies  from  0  to-. 23.   This  range  of  values 

r  D 
o  s 
of  —7 —  has  been  shown  to  have  only  a  minor  effect  on  the 

results. 

More  important,  as  shown  in  section  8,  it  is  basic  to 
the  tearing  mode  analysis  that  the  exponent,  h,  defined  by 
equation  (21a),  be  real.   Furthermore,  the  asymptotic 
analysis  requires  that  h  a-  0 .   Equation  (65)  shows  that  real 
h  is  guaranteed  for  6  <  .5,  while  Table  1  indicates  that 
reasonable  results  can  be  expected  from  the  asymptotic 
analysis  for  6  <  .2.   However,  we  carry  out  the  analysis  for 
values  of  5  as  high  as  .425  and  expect  the  results  to  have 
qualitative  validity. 

Exact  values  of  the  growth  rate  may  be  obtained  by 
using  the  following  procedure: 
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1.  Fix  the  position  of  the  outer  wall  at  some  point  at 
which  the  equilibrium  has  an  unstable  tearing  mode  for  the 
above  values  of  the  equilibrium  parameters  and  angular  velocity 
corresponding  to  6  =  .025.   As  mentioned  above,  the  criteria 
for  such  an  instability  is  that  A   >  0  where  A'  is  computed 
from  equation  (30b). 

2.  At  this  low  value  of  the  angular  velocity,  the  param- 
eter  A   governing  the  inner  equations  is  effectively  zero. 
Thus  the  criterion  of  merit,  F,  that  characterizes  the  effect 
of  Hall  current  and  rotation,  takes  a  value  approximately 
equal  to  its  vertex  value  of  F  =  2.12   as  shown  in  Figure  3. 
The  normalized  growth  rate,  y,    can  then  be  calculated  from: 

3/5  fAMl/5 


Y  =  n 


(*"] 


where  n  is  computed  from  (64). 

3.   The  angular  velocity  specified  by  6    is  increased  by 
a  small  amount  keeping  the  particle  number,  N  ,  fixed.   The 
outer  equations  are  reintegrated  to  obtain  a  new  value  of  A  . 
The  value  of  q  is  determined  from  equation  (62)  and  since  F 
depends  also  on  A   which  by  equation  is  a  function  of  the 
growth  rate  y,    the  growth  rate  must  be  determined  this  time 
from  equation  (61)  by  iteration.   This  iteration  is  initiated 
by  using  the  value  obtained  for  y    from  step  2  as  the  first 
guess . 
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4.    Step  3  is  repeated  for  the  case  of  no  Hall  current, 

i.e.,  e  =  q  =  0.   The  values  of  the  growth  rates,  are  used  as 

normalizing  factors;  i.e.,    normalized  =  ± 

Ye  =  0 

In  this  way  a  table  of  normalized  growth  rates  as  a 

function  of  angular  velocity  for  a  constant  wall  position 

would  be  obtained  for  one  family  of  equilibria  parametized 

by  NT .   There  is  however  one  problem  which  makes  the  im- 
J-j 

plementation  of  the  above  procedure  difficult.   The  deter- 
mination of  A'  requires  the  evaluation  of  (30b)  at  the  edge  of 
the  boundary  layer.   However,  as  shown  in  Table  1,  the  thick- 
ness of  the  boundary  layer,  x,  »  ,  as  calculated  by  equation 

(23),  varies  directly  with  angular  velocity  taking  on  values 
on  the  order  from  0  to  10   .   However,  this  close  to  the 
neutral  point,  where  B  (r  )  =0,  the  solutions  of  the  outer 
equations  are  asymptotically  large  and  the  integration  routines 
are  no  longer  sufficiently  accurate.   It  turns  out  that  for 
some  parameter  ranges  there  exists  a  small  annulus  about  the 
neutral  point  in  which  the  plasma  position  is  close  enough 
to  the  neutral  point  for  a  computation  of  i'  to  be  made 
yet  far  enough  away  from  the  neutral  point  to  enable  the 
outer  equations  to  be  accurately  integrated,  but  this  is 
generally  not  the  case. 

We  can  avoid  this  difficulty  of  integrating  the  outer 
equations  if  we  are  willing  to  settle  for  overestimates 
of  the  growth  rates  and  incomplete  information  as  to  the 
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exact  position  of  the  outer  wall.   We  merely  fix  the  value 
of  A   at  its  maximum  value  over  the  range  of  values  on  the 
interval  of  angular  velocities  of  interest.   This  corre- 
sponds to  some  position  of  the  outer  wall.   We  see  from 

equation  (61)  that  this  value  of  A   yields  an  overestimate 

-  -2 

of  y  whereas  from  equation  (63)  for  A  ,  overestimates  of 

-  -2 

Y  yield  smaller  values  of  A   for  the  same  value  of  q  thus 

from  Figure  3  this  results  in  lower  values  of  F  and  once 

again  from  equation  (61)/  overestimates  of  y-   We  have 

implemented  this  procedure  and  presented  the  results  in 

Figures  3~5,  where  ^-      is  plotted  against  the  normalized 

Ye  =  0 

angular  velocity,  6,  for  different  but  unknown,  outer  wall 
positions  characterized  by  and  varying  directly  with  A 
and  for  three  values  of  the  number  of  particles, 
NT  =  1016,  1017,  10   .   Actual  values  of  the  rigid  body 

Li 

angular  velocity,  fi,  and  ion  temperature  T  are  indicated 

in  figure  3. 

Figures  3~5,  show  that  the  inclusion  of  a  Hall  current 

results  in  increased  stability  to  tearing  modes  when  the 

plasma  column  rotates  in  one  direction  and  decreased  stability 

when  the  column  rotates  in  the  opposite  direction.   Figures 

3-5  show  further  that,  in  general,  both  stability  and  instability 

are  enhanced  by  rotating  the  column  faster  in  either  direction 

although  above  a  critical  value  of  the  wall  position,  A   , 

the  normalized  growth  rate,   —       , 

Ye  =  0 
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attains  a  maximum  with  respect  to  the  normalized  angular 
velocity,  5,  for  rotations  in  one  direction  and  a  minimum 
for  rotations  in  the  opposite  direction. 

Figures  3—5  also  show  that  the  growth  rates  are  very 
sensitive  to  the  wall  position,  A  ,  and  that  the  plasma 
column  becomes  more  stable  to  tearing  modes  for  outer  wall 
positions  nearer  to  the  plasma  vacuum  interface,  i.e.,  for 
smaller  values  of  A  .   Furthermore,  for  sufficiently  small 
values  of  A  ,  at  a  critical  value  of  the  normalized  angular 
velocity,  6,  tine  tearing  mode  instability  is  quenched,  i.e., 

X     =  0  . 

Ye  =  0 

Figures  3-5  also  show  that  for  a  given  wall  position, 
A  ,  the  stabilizing  effect  is  more  pronounced  at  lower  values 
of  the  number  of  particles,  N  ,  and  higher  ion  temperatures,  T, 

Li 
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7.   CONCLUSIONS 

Analysis  of  the  linearized  MHD  equations  of  motion  for 
a  rotating  field-reversed  theta  pinch  plasma  configuration  has 
shown  that  the  existence  of  Hall  currents  alters  the  dispersion 
realtion  of  unstable  tearing  modes  in  a  manner  which  is  consistent 
with  experimentally  observed  behaviour.   Specifically  the 
m=0,  k   7^  0  tearing  mode  growth  rate  can  be  reduced  for  preferred 
rotational  directions.   We  find  that  overestimates  of  the  effect 
of  Hall  current  and  rotation  on  the  growth  rate  of  tearing  modes 
can  be  determined  from  characterisitcs  of  the  equilibrium  at  the 
neutral  point  of  the  reversed  magnetic  field  provided  that  the 
equilibrium  is  unstable  to  tearing  modes ;   The  effect  of  the 

tearing  mode  is  characterized  by  the  criteria  of  merit   F   which 

—2 
depends,  tc  leading  order,  on  only  two  parameters   q  and  A 

When  q  <  0  the  Hall  current  always  enhances  the  stability  of  the 

plasma  system  to  tearing  modes  while  values  of  q  <  0  always  result 

_2 
in  reduced  stability.   The  parameter   A   is  governed  primarily 

by  the  position  of  the  outer  wall  and  the  value  of  the  angular 

velocity . 

For  a  given  number  of  particles,  the  stabilizing  effect  of  the 

Hall  current  is  enhanced  by  wall  positions  closer  to  the  plasma 

vacuum  interface.   For  wall  positions  closer  to  the  plasma  than 

a  critical  value,  the  tearing  mode  instability  is  quenched  at  a 

critical  value  of  the  angular  velocity  of  the  column.   For  wall 

positions  larger  than  this  critical  value,  the  growth  rate  attains 

a  minimum  at  one  value  of  the  angular  velocity  for  rotation  in 
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one  direction,  and  attains  a  maximum  for  another  value  of  angular 
velocity  for  rotation  in  the  opposite  direction.   For  a  given 
wall  position,  the  stabilizing  effect  of  the  Hall  current  is 
more  pronounced  for  fewer  total  particles,  NT  ,  and  higher  ion 
temperature,  T. 
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Figure  Captions 

Figure  la:   Analytical  solution,  £  (Z) ,  the  resistive  layer 

equations  for  a  special  case  (numerical  solution 
coincides  with  the  analytic  solution  within  the 
accuracy  of  the  graph) . 

Figure  lb:   Analytic  solution,  C(Z),  of  the  inner  equations  for 
a  special  case  (numerical  solution  coincides  with 
analytic  solution  within  the  accuracy  of  the  graph) . 

Figure  2:    Dependence  of  the  criteria  of  merit,  F,  on  the 

-2 

parameters  A   and  q  governing  the  resistive  layer. 

Figure  3:    Results  of  the  analysis  showing  the  effect  of  the 

normalized  growth  rate  —    on  the  normalized  angular 

Ye=0 

velocity  of  the  plasma  6,  and  a  measure  of  the  outer 

wall  position,  a',  for  particle  number,  N   =  10 

The  actual  angular  velocity  of  the  plasma  column  is 

noted  on  the  abscissa. 

Figure  4:    Results  of  the  analysis  showing  the  effect  of  the 

normalized  growth  rate  —    on  the  normalized  angular 

Ye  =  0 

velocity  of  the  plasma  5,  and  a  measure  of  the  outer 
wall  position,  A',  for  particle  number,  N   =  10 

Li 

The  actual  angular  velocity  of  the  plasma  column  is 
noted  on  the  abscissa. 
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Figure  5:  Results  of  the  analysis  showing  the  effect  of  the 

normalized  growth  rate  ^    on  the  normalized  angular 

Ye  =  0 

velocity  of  the  plasma  6,  and  a  measure  of  the  outer 

•  18 

wall  position,  A  ,  for  particle  number,  N   =  10 

Li 

The  actual  angular  velocity  of  the  plasma  column  is 
noted  on  the  abscissa. 
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